Resonant two-magnon Raman scattering in two-dimensional and ladder-type Mott 

insulators 
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We investigate the resonant two-magnon Raman scattering in the two-dimensional (2D) and 
ladder-type Mott insulators by using a half-filled Hubbard model in the strong coupling limit. By 
performing numerical diagonalization calculations for small clusters, we find that the model can 
reproduce the experimental features in the 2D that the Raman intensity is enhanced when the 
incoming photon energy is not near the absorption edge but well above it. In the ladder- type Mott 
insulators, the Raman intensity is found to resonate with absorption spectrum in contrast to the 2D 
system. The difference between 2D and the ladder systems is explained by taking into account the 
fact that the ground state in 2D is a spin-ordered state while that in ladder is a spin-gapped one. 



In Mott insulators, the nature of the photoexcited 
states is less clarified as compared with that of low-energy 
spin states that are described by the Heisenberg model. 
The excitations to the photoexcited states are observed 
by measurements of optical absorption. In addition, reso- 
nant Raman scattering is very useful for the investigation 
of the photoexcited states, because the resonance occurs 
when incoming photon energy uj\ is near the energy of 
the photoexcited states. 

In two-dimensional (2D) insulating cuprates that are 
typical Mott insulators, a two-magnon peak is observed 
in the Raman scattering |l| . Two-magnon is a magnetic 
excitation from the ground state of the Mott insulators. 
The iOi dependence of the two-magnon Raman intensity 
in the B lg geometry shows an interesting feature |^|, [|: 
The two-magnon Raman intensity is not enhanced at 
near the absorption edge (~ 1.8-2 eV), but a resonance 
occurs at higher u>\ (~ 3 eV). 

Two-magnon peaks are also observed in the ladder- 
type Mott insulators [Q. In contrast to the 2D sys- 
tem, the two-magnon Raman intensity in the ladder 
compound Sri4Cu2404i resonates with absorption spec- 
trum |BJ, i.e., the intensity is enhanced at the absorption 
edge (~ 2 eV). 

In this paper, we theoretically investigate the u>\ de- 
pendence of the two-magnon Raman intensity in the 2D 
and ladder-type Mott insulators, and we clarify the ori- 
gin of the difference of the u>i dependence between the 
2D and ladder systems. There have been several theoret- 
ical studies about the two-magnon Raman scattering in 
the ladder system, where only the non-resonant regions 
were examined || R E|. In contrast to them, our study 
deals with resonant regions. We find that numerical cal- 
culations for the Hubbard model in the strong coupling 
limit reproduce the difference of the w; dependence of 
the two-magnon Raman intensity at the absorption edge 
between the 2D and ladder- type Mott insulators. The 
difference is caused by the difference of spin states in the 
two insulators. 



The Hamiltonian of the Hubbard model reads 



ff=-*E (i 



Cj !<T + H.c. 



(1) 



where the summations run over nearest-neighbor 

pairs. We consider the strong coupling limit, i.e., U ^S> t. 
In this limit, the ground state \i) and the final states of 
the two-magnon Raman scattering |/) have no doubly 
occupied site, and they are described by the Heisenberg 
model with nearest-neighbor exchange interaction J(= 
At 2 /U). On the other hand, in the intermediate states 
of the Raman process, there are one doubly occupied site 
and one vacant site. In this case, an effective Hamiltonian 
of (fy) is given by 



U^UoHt^HtUo 



U. 



(2) 



where H t is the hopping term of H, and Ho, III, and H2 
are projection operators onto the Hubert space with zero, 
one, two doubly occupied sites, respectively. A complete 
expression of Eq. (||) has been given elsewhere . The 
value of U is evaluated to be U=lOt for the gap values 
to be consistent with experimental ones. 

The Raman scattering intensity at the Ra- 
man shift energy u> is expressed as R(u>) — 
J2 f \(f\M R \i)\ 2 5{Lu-E f +E z ), where E t and E f 
are the ground-state and final-state energies, respec- 
tively. Mr is the Raman tensor operator whose matrix 
element is given by 



(i|^l/> = E 



{f\j a \y){y\jfi\i) 
En— Ei- tot + iT ' 



(3) 



where a and /3 are polarizations of the incoming and 
scattering photons, respectively, j a is the current oper- 
ator to the a direction and, E^ and T are energy and 
a damping rate for the intermediate states, respectively. 
Absorption spectrum £2 for the a polarization is given by 
^ 2 E M K*|ia|M>l a *(^ - E„ + Ei). We calculate R(u) 
and £2 by using the numerical diagonalization method 
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FIG. 1: The u>i dependence of Raman intensity i?2ma g (dashed 
hne), and absorption spectrum £2 (solid line) in a 20-site 2D 
cluster with U = lOt. Inset shows two-magnon Raman spec- 
trum R (ui) at uii — 4.0t . In £2 and R (lu), the solid lines are 
obtained by performing a Lorentzian broadening with of QAt 
on delta functions denoted by thin vertical bars. 



for small clusters (V20 x v20 square and 10 x 2 ladder 
lattices) . We use periodic boundary conditions in the 2D 
system. In the ladder system, a periodic boundary con- 
dition is used along the leg direction (the x direction), 
while an open boundary condition is used along the rung 
direction (the y direction). 

First, we show results in the Bi ff geometry for the 2D 
system in Fig. [j] |Tc| ]. The inset shows R (to) at lo\ — 4.0t. 
We define i?2mag as the integrated intensity of the peaks 
at < lu < 2.0t. The u>i dependence of i?2ma g (dashed 
line) and £2 (solid line) are shown in the main panel. £2 
has an edge peak at u>i = 6.2t. However, i?2mag is not 
enhanced at the peak. This lu- x dependence of i?2mag is 
similar to the experimental data |j| . 

Next, we discuss results for the ladder-type Mott in- 
sulators. We calculate R(lu) for the (xx) polarization, 
where (xx) means that the polarizations of both the in- 
coming and scattering photons are along the x direction. 
The values of hoppings (exchange interactions) are taken 



to be 



= t 



rung 



(Jieg = Jrung)- We note that the 



following results are qualitatively unchanged even if we 
use a realistic rate, J rU ng/ ^lcg = 0.8, corresponding to 
Sri 4 Cu 2 404i [||. In Fig.g, the inset shows R(lu) for the 
(xx) polarization at LUi—A.0t. The spectrum spreads in 
the range of 0.5t < lu < 1.5t, but the main peak ap- 
pears at lu ~ l.Oi. This is consistent with the fact that 
the peak obtained by the non-resonant Raman scattering 
calculations Jf], || is located at lu ~ l.Ot. The depen- 
dence of i?2mag (dashed line) for the (xx) polarization 
and £2 (solid line) for the x polarization is shown in the 
main panel. £2 has an edge peak at Wi=6.5t. i?2mag is 
enhanced at the peak in contrast to the 2D system, being 
consistent with experimental features. 

Let us consider the origin of the different Wi depen- 
dence between the 2D and ladder systems. In the 2D 
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FIG. 2: The u>i dependence of Raman intensity i?2ma g (dashed 
line), and absorption spectrum £2 (solid line) in a 10 x 2 ladder 
cluster with U = lQt for the (xx) polarization. Inset shows 
two-magnon Raman spectrum 7? (lu) at lu[ = 4.0t. In £2 and 
R (lu), the solid lines are obtained by performing a Lorentzian 
broadening with of OAt on delta functions denoted by thin 
vertical bars. 



Mott insulators, the ground state \i) is dominated by 
Neel-type spin configurations because of the presence of 
antiferromagnetic order, while the weight of the Neel- 
type spin configurations in the two-magnon final state 
I/) is very small JToj ] . Since the photoexcited state at 
the absorption edge, \[i dgc), significantly contains the 
Neel-type spin configurations in the spin background, 
(Mcdgoj jx \i) becomes large but (f\j x |^ cd gc) is very small 
because of a small amount of the Neel-type spin config- 
urations in I/). Therefore, i?2mag, which is proportional 

to K/l 3x l^cdgc) (Mcdgd jx , is not enhanced at the ab- 
sorption edge. On the other hand, in the ladder system 
we find that the weight of the Neel-type spin configura- 
tions in \i) is smaller than that of 2D but the weight in 
I/) at lu — 0.6t is as large as that in |z) and is thus much 
larger than that of 2D. As a result, (f\j x l/^edge) is en- 
hanced as compared with that in 2D and thus i?2mag is 
enhanced at the absorption edge. This is intuitively un- 
derstood in the following way. The two-leg ladder has a 
spin gap and thus the ground state is a spin-liquid state in 
contrast with the 2D system where antiferromagnetic or- 
der exists. While in 2D the two-magnon final state is very 
different from the initial state with respect to the Neel- 
type spin spin configurations, the difference between the 
initial and final states in ladder is small because both the 
states are of spin liquid and contain a similar amount of 
the Neel-type spin configurations. This clearly explains 
the different lu\ dependence between 2D and ladder. 

Finally, we show results of the two-magnon Raman 
spectrum for the (yy) polarization, i.e., the poralization 
along the rung direction. In Fig. [|, the inset shows R (lu) 
at LUi — LOt. The spectral shape is almost the same as 
that for the (xx) polarization at lu\ = 4. Of. We note that 
in the calculation of the non-resonance Raman scattering 
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FIG. 3: The same as Fig. q but for the (yy) polarization. 

the same shapes are obtained for the (xx) and (yy) po- 
larizations in the case of J\ cg — J rung . The main panel 
in Fig. U shows the u>i dependence of i?2mag (dashed line) 
for the (yy) polarization and £2 (solid line) for the y po- 
larization. £2 exhibits an edge peak at uj\ — 7.3t. i?2mag 
is resonantly enhanced at oj\ ~ lOi, but not at the ab- 
sorption edge in contrast to the (xx) polarization. This 
is, however, in contradiction with the experiment where 
#2niag f° r (yy) is enhanced at the absorption edge [gj. 
This disagreement may come from the fact that we do 
not consider the interactions between ladders from the 
theoretical viewpoint and a small amount of holes are 
doped into the ladders in Sri4Cu2404i from the experi- 
mental side. 

In summary, we have investigated the dependence of 
the two-magnon Raman intensity on the incoming pho- 
ton energy uj[ in the 2D and ladder-type Mott insula- 



tors by using the Hubbard model in the strong coupling 
limit. We have found that the difference of spin states 
between the 2D and ladder systems causes different LOi 
dependence of the two-magnon Raman intensity at the 
absorption edges. 
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